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An  algorithm  is  presented  for  solving  minimum- cost  flow 
problems  in  which  each  arc  of  the*  network  bus  a  finite* 
maximum  flow  capacity  and  a  concave  cost  function  associated 
with  sending  flow  alonq  that  arc.  each  cost,  function  is 
broken  into  a  series  of  cost  increments  throuqh  the  use*  of 
piecewise  linear  approximat ions.  The  algorithm  taken  any 
feasible  solution  and  recirculates  flow  o\"-r  loss  costly 
cycles  to  obtain  an  optimal  solution.  A  pod  if ication  which 
handles  the  existence  of  non-zero  lower  bounds  on  flow 
through  the  various  arcs  is  also  yiven. 


/ 


DD 


'.V.-..1473 


(PAGE  1) 


UKCLAShli  T!  0 


S/H  O’  01  •«07*MI  I 


Uiunij  CU«»ilicalt>>n 


Mlltl 


Network  Flown 

Optimization 

Concavity 


it 


DD  ,"'.‘..1473  <e*cK. 


Minimum-Cost  Flown  in  Networks 
with 

Upper  Doundod  Area  and  Concave  Coat  Functions 


Author 


Approved 


by 


Wayne  Jay  Hnllonbock,  Jr. 
Lieutenant  Junior  Grade,  United  Statoa  Navy 
D.S.,  United  States  Naval  Academy,  1971 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  1972 


sL'-.'L  j  .  X\ 


by: 


~  Ld  /fhYi,  J7f _ , _ 

Thesis  Advisor 


Second  Reader 


_ _ ee  /  <  '  _ 

"cRalrman,  Department  GT  Operations  Research 
J  and  Administrative  Sciences 

*/.•»'  .*  •  t  / ^  , 

_ 

‘  71  '  Academic  Dean 


/ 


TABLE  OF  CONTENTS 

I.  INTRODUCTION . - . . . .  6 

A.  TAST  WORK - 6 

D.  OBJECTIVE  AND  SCOPE - - - - -  7 

II.  MODEL  DEVELOPMENT - - -  9 

A.  NETWORK  DESCRIPTION.*:  AND  ASSUMPTIONS -  9 

B.  FORMULATION  OF  THE  COST  FUNCTION - -  10 

C.  STATEMENT  OF  THE  COST  MINIMIZATION  PROBLEM--  J1 

D.  LENGTH  SEQUENCE  NETWORKS -  11 

E.  DEFINITIONS  OF  CYCLES -  15 

III.  SOLUTION  PROCEDURE -  17 

A.  PREVIEW - -  17 

B.  THE  ALGORITHM -  18 

IV.  EXAMPLE  PROBLEM -  20 

V.  EXTENSIONS - - 25 

A.  SPECIAL  CASE:  Q  1  MINIMUM  f$i  .  VALUE -  25 

B.  NON-ZERO  LOWER  BOUNDS -  20 

VI.  AREAS  FOR  FURTHER  S'lUDY -  28 

A.  IDENTIFICATION  OF  NEGATIVE  CYCLES -  28 

B.  GENERAL  COST  FUNCTIONS -  20 

VII.  SUMMARY - - 30 

BIBLIOGRAPHY -  31 

INITIAL  DISTRIBUTION  LIST - - -  32 

FORM  DD  1473 . . - - - - -  33 


3 


LIST  or  ILLUSTRATIONS 

Figure 

1.  Construction  of  a  Length  Sequence  Network  from 


a  Feasible  Solution  of  a  Given  Network  Problem——  14 

2.  Examples  of  Cycles - 16 

3.  Example  Network - 20 

4.  Initial  Feasible  Solution — - - 20 

5.  Initial  Length  Sequence  Network - 21 

6.  First  Reallocation  of  Flow - 22 

7.  Second  Length  Sequence  Network - - 22 

8.  Second  Reallocation  of  Flow - 23 

9.  Third  Length  Sequence  Network - - -  24 


1 


I 

! 

1 


4 


tabu:  OF  SYMBOLS  AMI)  AUbKLVIATIONS 


G  (N«A) 

(i#  j) 

Mij 

Lij 

xii 

cij<v 

C(K) 

cij 


0 

z 

X 


dii 

y(k) 


V 
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a  connected  original  network  consisting  of  a 
sot  of  nodes  N  and  a  nui  of  arcs  A 

an  arc  directed  from  node  i  to  nodo  j 

the  upper  bound  on  the  flow  capacity  of  arc 
(i,j)  of  the  original  network 

the  lower  bound  on  the  flow  capacity  of  arc 
(i,j)  of  the  original  network 

thu  actual  amount  (an  integer  value)  of  flow 
in  arc  (i,j)  of  the  original  network 

the  cont  of  sending  X^  unit:*  of  flow  along 
arc  (i,j)  of  the  original  network 

the  incremental  cost  of  tending  a  kth  unit  of 
flow  along  arc  (i,j)  of  the  original  network, 
given  that  k-1  unite  have  already  been  sent 

the  amount  of  flow  to  bo  sent  through  the 
original  network 

the  total  cent  of  sending  Q  units  of  flow 
through  the  original  network 

a  flow  vector  having  all  of  the  various  X^. 
values  as  its  components  J 

a  connected  length  sequence  network  consisting 
of  a  set  of  nodes  N  and  a  set  of  arcs  A| 

the  kth  element  in  the  length  eccjucnco  of  arc 
(i,j)  of  tho  length  sequence  network 

tho  length  of  the  directed  arc  from  node  i  to 
node  j  in  o  shortest  route  algorithm 

tho  label  value  of  node  i  upon  completion  of 
tho  kth  iteration  of  Dijkstra's  algorithm 


tho  permanent  label  value  of  nodo  i 


1 •  INTRODUCTION 

A.  PAST  WORK 

Considerable  literature  exists  dealing  with  minimum-cost 
network  problems  in  which  the  various  arcs  of  a  network  have 
upper  and  lower  bounds  on  flow  capacity  and  1  inoar  cost 
functions.  In  particular,  thoro  is  the  Primal-Dual  Algorithm, 
developed  by  Ford  and  Fulkerson  [2]  in  1955,  which  solves  such 
network  flow  problems  having  all  lowor  bounds  equal  to  zero. 
Ford  and  Fulkerson  (5)  later  developed  the  Out-of-Kilter 
Algorittim  to  solve  problems  having  non-zero  lower  bounds. 

The  two  above  algorithms  can  also  be  used  to  deal  with 
minimum-coot  network  problt?ma  in  which  the  cost  functions 
are  convex.  By  making  piecewise  linear  approximations  to 
the  convex  cost  curves,  every  arc  in  such  a  network  can  be 
replaced  by  a  group  of  arcs  having  different  linear  cost 
functions,  where  each  arc  of  the  group  corresponds  to  one 
sogmont  of  tho  linear  approximation  to  the  convex  cost 
function. 

liu  (6)  has  looked  at  a  special  case  of  using  piecewise 
linear  approximations  in  convex-cost  networks.  Ho  defines 
"up"  and  "down"  arc  costs  os  tho  incremental  costs  incurred 
from  increasing  or  decreasing,  respectively,  by  one  unit  of 
flow,  the  already  oxisting  flow  in  an  arc.  He  presents  a 
solution  procedure  which  optimally  increments,  from  zero  to 
any  desired  amount,  the  total  flow  through  the  network. 
Although  Hu  does  not  consider  upper  )>ounds  on  arc  flow 
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capacity,  hiu  alyorithm  can  easily  he  modified  to  bundle 
thin  condition  by  making  the  "up-eont"  of  a  Maturated  arc 
infinite. 

Itowovor,  in  many  practical  application!;  the  convex  coot 
assumption  does  not  hold,  bant/ i ij  (1)  pointu  out  that  the 
presence  of  a  oot-up  or  red  tape  charyo  yields  a  concave 
coot  function.  Likewise,  it  can  bo  sryued  that  efficiencies 
of  scale  and  tho  practice  of  yiviny  discounts  or  rebates 
in  transact  ions  involving  larye  quantities  of  cjoodu  or  ser- 
vicoo  also  yield  concave  coot  functions.  In  yonoral,  it  can 
bo  stated  that  concave  co«t  functions  arise  through  tho 
existence  of  decreasing  naryinal  costs,  a  common  phenomenon 
in  roal  life  situations.  Thus,  there  is  yood  reason  to 
study  concave  cost  networks. 

Zangwill  (0)  presents  a  solution  to  the  minimum-cost 
flow  problem  in  which  tho  various  arcs  of  a  network  have 
concave  cost  functions,  althouyh  he  does  not  consider  the 
existence  of  upper  or  non- aero  lov/er  bounds  on  the  arc 
capacities.  Tho  solution  to  his  problem  has  all  of  tho 
flow  bciny  sent  alony  tho  minimum  total  cost  chain  of  arcs 
from  tho  source  to  the  sink. 

B.  OBJECTIVE  AND  SCOPE 

The  pur  peso  of  this  paper  is  to  present  a  method  of 
solviny  minimum-cost  flow  problems  in  which  tho  arcs  of 
tho  network  have  zero  lower  bounds,  finite  upper  bounds, 
and  concave  cost  functions. 


7 


For  any  given  amount  of  flow  which  must  bo  sent  from 
the  source  to  tho  sink,  tho  problem  can  be  solved  if  the 
maximum  flow  capacity  and  the  cost  function  are  known  for 
each  arc  in  the  network.  A  basic  assumption  to  the  solution 
procedure  is  that  oach  cost  function  can  be  broken  into  a 
series  of  non-increasing  cost  increments  through  the  use  of 
piecewiso  linear  approximations. 

As  outputs  to  the  problem,  the  solution  procedure  pro¬ 
vides  tho  optimal  arc  flows  and  the  associated  optimal 
total  cost. 


II.  MODEL  DEVELOPMENT 


A.  NETWORK  DESCRIPTION  AND  ASSUMPTIONS 

Consider  a  connected  network  G (N , A)  consisting  of  a 
set  of  nodes  N  and  a  set  of  arcs  A.  Let  the  integers 
i=l,2,,..,n  represent  the  nodes  and  the  two-tuples  (i,j) 
(i=l, 2, . . . ,n;  j=l,2,...,n;  and  i  ^  j)  represent  the  arcs. 
Let  node  1  correspond  to  the  source  and  node  n  correspond 
to  the  sink.  The  arcs  are  assumed  to  be  directed  so  that 
the  order  (i,j)  implies  an  arc  directed  from  node  i  to 
node  j . 

It  is  assumed  that  there  is  only  one  source  node  and 
only  one  sink  node.  It  is  also  assumed  that  no  more  than 
one  arc  connects  any  pair  of  nodes  in  the  same  ordered 
direction.  If  multiple  sources  or  sinks  exist,  artificial 
nodes  and  arcs  may  be  added  to  the  network  so  that  there 
is  a  single  overall  source  and  a  single  overall  sink. 
Similarly,  if  two  or  more  directed  arcs  connect  the  same 
pair  of  nodes  in  the  same  ordered  direction,  artificial 
nodes  and  arcs  may  be  added  to  the  network  as  needed  tc 
alleviate  this  condition. 

Associated  with  each  arc  are  an  upper  bound  on  capa¬ 
city  and  a  cost  function  for  flow  over  the  arc.  Let 
be  the  actual  amount  of  flow  and  let  M. .  be  the  maximum 
flow  capacity  for  an  arc  (i,j).  Then,  for  to  be  a 
feasible  arc  flow, 
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for  each  arc  {i,j)  in  the  network. 

Conservation  of  flow  is  assumed  to  exist  at  all  nodes. 
The  total  flow  out  of  any  node  equals  the  total  flow  into 
that  node;  that  is,  there  is  no  storage  of  flow  at  any  of 
the  nodes. 


B.  FORMULATION  OF  THE  COST  FUNCTION 


Let  the  concave  cost  function  CjLj(X^j)  represent  the 


cost  of  sending  an  amount,  X^ j ,  of  flow  along  arc  (i,j). 


It  is  assumed  that  C-j-j(Xj.)  is  continuous  and  non-negative 
over  the  entire  range  (from  0  to  M^j)  of  X^j.  It  is  also 


assumed  that  C^ ^  (0)  •=  0  for  each  arc  in  the  network.  How¬ 


ever,  if  there  is  a  cost  function  such  that  C^ j  (0)  ^  0, 


then  a  new  cost  function  C^j  (X  — )  =  C^jfX^j)  -  C^j  (O)  , 


can  be  substituted  for  C^j  (X-^j)  without  changing  the 
nature  of  the  optimization  problem. 

It  is  assumed  that  piecewise  linear  approximations  can 
be  made  to  the  cost  function  of  every  arc  in  the  network. 

If  the  arc  flows  and  arc  flow  capacities  are  required  to 
be  non-negative  integers,  the  cost  function  of  an  arc  (i,j) 
can  be  replaced  with  a  non-increasing  sequence  of  incremental 
costs;  C.; ,  C_.  j  , .  . . ,  Cj :  ,  where 


(1) 

(2) 

c  (Mi5 
•  * '  Lij 

cij 

'  Cij  '* 

c!H)  =  ci}W 

-  (k-l) 

(2) 


for  k  =  1,2, .. . /M^j . 
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Thu  interpretation  hero  iu  that  an  incremental  cost, 

,  in  tho  additional  coot  of  sending  one  unit  of  flow 
across  arc  (i,j),  given  that  k-1  unitu  of  flow  have  already 
boon  sent  across  arc  (i,j).  Thus,  if  X^j  represents  the 
actual  amount  of  flow  sent  across  an  arc  (i,j),  it  can 


oasily  bo  seen  that  the  associated  cost, 

(0)  (X,,1 

Cij  *xij*  "  Ci  j  +  ‘  *  ■*'  cij  • 


whore  is  defined  to  be  equal  to  zoro. 


C.  STATEMENT  OF  THE  COST  MINIMIZATION  PRORLEM 


All  flow  through  the?  network  in  assumed  to  travel  froii 
the  source  node  to  the  sink  node.  Then,  for  any  specified 
amount  of  flow  Q  to  be  sent  through  the  network,  the  cost 
minimization  problem  can  be  stated  as  the  following  pro¬ 
gramming  problem: 


Find  non-negative  integer  values  for  all  X^  which 


minimize  Z  =  £2 

U,j>  €  A 


( s  •»")  ■ 


subject  to  Tj  X,  .  -  Zj  X. 


Q,  i*-l 

0,  i=2,3, . . . ,n-l 

-0,  i=n  ^5) 


and  0  <  Xi;.  <  Mi_. ,  VU.jUA.  (6) 

D.  LENGTH  SEQUENCE  NETWORKS 

Let  X  denote  a  flow  vector  which  has  as  its  components 
the  X..  values  of  all  the  arcs  in  the  network.  Any  flow 
vector  X  which  satisfies  constraints  (5)  and  (6)  of  the 
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coot  minimization  problem  in  a  fo.oiihlc  solution  to  the  pro¬ 
blem.  The  net  o(  all  (low  vac  torn  which  are  feasible 
solutions  comprise  the  problem'll  feasible  region. 

Consider  any  flow  vector  X  which  it  a  feasible  solution 
to  Lho  problem.  Associated  with  thin  solution,  a  unique 
connected  length  sequence  network  C|  (N,A|  )  can  be  con¬ 
structed.  This  new  network  has  the  same  nodes  as  the 
original  network,  but  has  a  different  arrangement  of  arcs, 
as  described  bclowt 

1. )  If  arc  (i,j)  of  the  original  network  is  empty,  draw 
a  forward  arc  from  node  i  to  node  j  in  the  length  sequence 
network.  Label  this  forward  arc  with  the  length  sequence 

(1)  CM.  )  (k)  <k) 

ij  . 'ij  3  -  whuro  »  c4j  . 

2. )  If  arc  ( i , j )  of  the  original  network  is  neither 
empty  nor  saturated,  chaw  a  forward  arc  fro».  node  i  to  node 
j  and  a  backward  arc  from  node  j  to  node  i  in  the  length 
sequence  network.  Label  the  forward  arc  with  the  length 


coqucnco 


I 


(k) 

ij 


C 


<vk) 

ij 


•  •  I 


I 


("ij-Xij 

ij 


,  where 


Label  the  backward  arc  with  the  length 


CD  l(xi1)  .00  CX^  +  l-k) 

soquoncc  »•••*  *ji  J  »  where  ■  -C^  J 

3.)  If  arc  (i,j)  of  the  original  network  is  saturated, 
draw  a  backward  arc  from  node  j  to  node  i  in  the  length 
sequence  network.  Label  this  backward  arc  with  the 
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length  soquoneo 


(1)  (M,,) 

*ji . Iji  J 


who  i  t? 


<k)  (M,,+l-k) 

»n  *VJ 


Lot  tho  pair  of  "forward"  and  "backward"  arc*  corres¬ 
ponding  to  thu  uomo  original  network  arc  bo  called  a  com** 
plcmont .  Tho  number  of  complements  in  a  length  sequence 
notwork  is  equal  to  tlio  number  of  area  in  tho  corresponding 
original  network  which  arc  noitln  r  empty  nor  saturated. 

Figure  1  illustrates  how  a  length  sequence  network  is 
constructed  from  and  corresponds  to  i  feasible  solution  of 
a  simple  notwork  flow  problem  having  0  ■■  3.  Note  tho  two 
complements  £(1,2) ,  (2,1)J  and  |  (1,3) , (3, 1)J  in  the 
length  sequence  notwork  corresponding  to  tho  two  arcs  (1,2) 
and  (1,3)  in  the  original  network  which  are  neither  empty 
nor  saturated. 

For  a  forward  arc  (i,j)  in  the  length  sequence  network, 
tho  sequence  of  length  numbers  corresponds  to  the  sequence 
of  incremental  costs  incurred  with  sending  additional  units 
of  flow,  in  excess  of  X^j,  along  the  unrnturatod  arc  (i,j) 
of  the  original  network.  Similarly,  for  a  backward  arc 
(j,i),  the  sequence  of  longth  numbers  corresponds  to  tho 
sequence  of  incremental  costs  incurred  with  "unsending,"  or 
removing,  existing  units  of  flow  from  the  non-empty  arc 
(i,j)  of  the  original  network.  Duo  to  the  concave  nature  of 
the  cost  functions  involved,  tho  sequence  of  length  numbers 
associated  with  each  arc  in  the  length  sequence  network  is 
non- increasing. 
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Figure  1.  Conutruction  of  n  Length  Sequence  Network  fron 
a  Feasible  Solution  of  a  C’.ivun  Network  Problem 


(1)  (2)  (3) 


(c)  Correnponding 
Length  Sequence 
Network 


k .  Di-:piNiTiot:s  or  cyci.i::; 

With  respect  to  a  length  sequence  network,  a  simple 
cyclo  will  bo  defined  as  a  series  of  adjacent  and  consist- 
ontly  directed  area  which: 

(1)  begins  and  onda  at  tho  same  noJo, 

(2)  passes  through  at  least  two  other  nodes, 

(3)  doou  not  havo  any  arc  in  the  sorion  which  itt 
traversed  more  than  once,  and 

(4)  doou  not  include  both  member b  of  any  complement. 

If  each  arc  of  a  Bionic  cycle  in  a  length  sequence  net** 

work  hue  at  leant  two  olementu  in  its  length  sequence,  then 
it  is  possible  to  expand  the  simple  cycle  into  a  double 
cycle  by  identically  and  complete  ly  repeating  tho  sinplo 
cycle.  Similarly,  if  each  arc  of  a  simple  cyclo  in  a 
length  sequence  network  has  at  least  three  elements  in  its 
length  sequence,  then  it  is  possible  to  expand  the  simple 
cycle  into  a  triple  cycle  by  identically  and  completely 
repeating  tho  uinplc  cycle  twice.  A  multiple  cycle  will 
bo  defined  as  any  double  cyclo,  triple  cycle,  ot  cetera. 

A  compound  cycle  will  be  defined  to  be  a  combination 
of  two  or  more  simple/nultiplo  cycles  such  that: 

(1)  each  simplc/multiplc  cycle  in  the  combination  has 
at  loast  ono  arc  in  common  with  another  simplc/multiplo 
cycle  in  the  combination, 

(2)  each  arc  common  to  two  or  more  simplc/multiplc 
cycles  has  a  sufficient  number  of  elements  in  its  length 
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•  eotjuencu  to  accommodate  dll  of  the  .iiinplc/»nul  t.  ipic  cycles 
to  which  it  in  common,  and 

(3)  not  more  than  one  ai/c  of  any  complement  in  in¬ 
cluded  in  the  compound  cycle. 

The  term  cycle  will  hereafter  bo  uuod  to  roprenent  any 
of  the  specific  cyclon  defined  above.  Figure  2  illustrates 
examples  of  cyclos. 
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(1)  (2) 
21  '*21 


Simple  Cycles: 
Doublo  Cycle: 
Triplo  Cycle: 
Compound  Cyclos: 


1-3-4-2-1;  1 - 3 - 2 - 1 ;  2-3-4-2 

1- 3-4-2-1- 3-4-2-1 
None  exist 

1  —  3— 4— 2— 1  —  3— 2—1;  2- 1 - 3-4 - 3-4 - 2 ; 

1-3-4-2-1-3-4- 2-1- 3-2-1 


Noto  that  2-3-4-2-1-3-4-2-1-3-2  is  not  a  legitimate 
cyclo  since  it  uses  both  memberu  of  tho  complement 
{(2,3). (3,2)}. 
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III.  HOI.U1- 1 0!;  I '  Hi: 


a.  i*iu:viki; 

Tho  algorithm  to  bt>  prosonud  begins  by  dcl.-i ui ni ng  any 
feasible  solution  to  tho  given  network  problem.  A  length 
sequence  network  can  then  be  construct « d  as  described  above, 
Tho  identification  of  cycles  in  the  length  sequence  network 
which  havo  negative  total  length  results  in  a  ^circulation 
of  flow  in  the  original  network  at  a  reduced  total  cout. 

Thin  recirculation  ol  flow  results  in  a  char.y*  of  the  exiut- 
ing  feasible  solution.  A  new  length  sequence  is  then  created. 
Tho  identification  of  cycles  having  negative  length  in  thin 
new  length  sequence  network  trigger.*!  another  iteration.  The 
algoritlim  terminates  when  no  cycles  having  negative  length 
can  bo  identified  in  the  length  sequence  netv/ork  corres¬ 
ponding  to  the  existing  feasible  solution  of  the  original 
network.  The  existing  feasible  solution  at  the  time  of 
termination  is  an  optimal  solution  to  the  problem. 

The  algorithm  itself  docs  not  identify  cycles  in  the 
length  sequence  network  which  have  negative  total  lengths. 

Due  to  tho  largo  number  of  nodes  and  arcs  which  may  exist 
in  a  length  sequence  network,  identification  of  such  cycles 
may  be  a  difficult  task.  In  order  to  terminate  the 
algorithm,  all  possible  cycles  in  the  length  sequence  net¬ 
work  must  be  chocked  to  ensure  that  the  existing  feasible 
solution  to  the  problem  is  an  optimal  solution. 
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H .  TIIK  AI.COKITHM 


X.  Kim',  any  feasible  solution  to  lliu  problem  and  send 
flow  acrous  the  network  accordingly.  Use  (4)  to  determine* 
the  total  trout  associated  with  thin  initial  feasible 
solution. 

2.  based  on  the  existing  feasible  flow,  confitruct  a 
length  sequence  network  using  tho  steps  described  in  Section 
III). 

3.  Try  to  find  a  cycle  in  tho  length  sequence  notwor); 
which  has  a  negative  length.  In  attempting  to  identify 
such  a  cycle,  two  important  conditions  must  bo  met: 

a.  If  a  particular  arc  in  tin?  length  sequence 
network  is  traversed  k  times,  then  the  first  k  elements  in 
that  arc's  length  sequence  must  be  used  in  computing  the 
cycle's  total  length. 

b.  Not  more  than  one  arc  of  any  complement  nay  be 
traversed. 

4a.  If  a  cycle  of  negative  length  can  be  identified, 
reallocate  flow  in  tho  original  network  as  follows: 

(1)  If  forward  arc  (i,j)  of  the  length  sequence 
network  is  traversed  k  times  by  the  cycle  of  negative 
length,  increase  the  flow  through  arc  (i,j)  of  the 
original  network  by  k  units. 

(2)  If  backward  arc  (j,i)  of  the  length  sequence 
network  is  traversed  k  times  by  the  cycle  of  negative 
length,  decrease  the  flow  through  arc  (i,j)  of  the 
original  network  by  k  units. 
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Also,  dotcimino  tho  total  coot  atmociatud  with  this 
now  feasible  flow.  This  will  bo  equal  to  the  caul  associ- 
otod  with  the  previous  feasible  flow  plan  the  length  of  tho 
cyclo  which  hats  just  been  identified. 

Return  to  step  2. 

4b.  If  a  cycle  of  nogAtivo  length  cannot  be  found,  then 
terminate  tho  algorithm.  The  existing  foasible  flow,  so  well 
as  tho  associated  total  cost,  is  optimal.  If  thoro  oxiots 
a  cyclo  of  length  zero,  thcro  is  an  alternate  optimal  solu¬ 
tion  to  tho  problem  obtainable  by  carrying  out  step  4a. 
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iv.  i:xahm:  phohi.*  m 

Consider  Hit*  network  shown  in  Figure  ,1.  *i)*o  numerical 

quantities  associated  with  each  arc  (i#j)  arc: 

(1)  (2)  (M,,) 

Mij  (Cij  »Cij  *  * 


The  alogrithn  presented  above  will  bo  uretl  to  solve 
this  problem  for  a  specified  Q  ■  4. 

Stop  1.)  Fiyurc  4  shows  an  arbitrary  initial  feasible 
solution.  The  total  cost  associated  with  this  initial 
feasible  solution  is:  Z  -  $55. 


20 


Initial 

Feasible 

Solution 


Step  2.)  Figure  %  allows  the  length  sequenco  network 
corresponding  to  this  initial  feasible  flow. 


Figure  5.  Initial  Length  Sequence  network. 

Step  3.)  Tho  simple  cycle  2-5-6“4-3~l  -2  ban  a  total 

(!)  U)  (1)  (1)  (1)  #(1) 

length  equal  to  1 ^  ♦  1 «  {43  ♦  "31  4  *\2  m  4C* 

"Vo  complete  cyclcH  via  2-S-6-4- 3-1-2  have  a  length  of 
♦6  ♦  0  •  ♦ 6 .  Three  complete  cycles  via  2-5-6-4-3-1-2  have 
a  length  of  +6  ♦  0  -8  «  -2.  Therefore,  a  cycle  of  negative 
length  huu  been  found. 

Step  4a.)  Forward  arcs  (1,2), (2,5),  and  (5,0  of  tho 
length  sequence  network  are  each  traversed  three  limes  by 
the  cyclo  of  negative  length  which  has  just  boon  identi¬ 
fied.  Thorofore,  tho  flow  through  nren  (1,2),  (2,5),  and 
(5,6)  in  tho  original  network  is  increased  by  three  units, 
backward  arcs  (6,4),  (4,3),  and  (3,1)  of  tho  length  sequence 
network  arc  each  traversed  three  times  by  the  cyclo  of 
negative  length  which  has  boon  identified.  Therefore,  tho 
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flow  through  arcs  (4,C),  (3,4),  mill  (1,3)  in  tho  original 
notwork  iu  decreased  by  throe  unite.  Figure  6  ohowe  tho 
resulting  real  location  of  flew. 


3 


The  total  cost  associated  with  this  new  feasible 
solution  is:  Z  ■  55  -  2  ■  $D3. 

Stop  2.)  Figure  7  shows  the  length  sequence  network 
corresponding  to  this  new  feasible  flow. 


Figure  7.  Second  I.ength  Sequence  Network. 
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Step  3.)  The  simple  cycle  2-5-4-3-2  in  this  second 
length  sequence  network  has  a  total  length  equal  to 

(1)  (1)  (1)  (1) 

J  25  +  * 54  +  4  43  +  *32  =  “7 *  Therefore,  a  cycle  of  nega¬ 

tive  length  has  been  found  in  this  second  length  sequence 
network.  Unfortunately,  this  cycle  cannot  be  repeated. 

Step  4a.)  Since  the  forward  arc  (2,5)  of  the  length 
sequence  network  is  traversed  one  time  by  the  recently  dis¬ 
covered  cycle  of  negative  length,  the  flow  through  arc  (2,5) 
in  the  original  network  is  increased  by  one  unit.  Similarly, 
since  the  backward  arcs  (5,4),  (4,3),  and  (3,2)  of  the  length 
sequence  network  are  each  traversed  one  time  by  the  recently 
discovered  cycle  of  negative  length,  the  flow  through  arcs 
(4,5),  (3,4),  and  (2,3)  in  the  original  network  is  decreased 
by  one  unit.  Figure  8  shows  the  resulting  reallocation  of 
flow. 

4 


0 
5 

Figure  8.  Second  Reallocation  of  Flow. 

The  total  cost  associated  with  this  new  feasible 
solution  is:  Z  =  53  -  7  =  $46. 
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Step  2.)  Figure  9  shows  the  length  sequence  network 


corresponding  to  this  nev;  feasible  flow. 


Figure  9.  Third  Length  Sequence  Network. 

Step  3.)  There  are  no  cycles  of  negative  length  in 
this  third  length  sequence  network. 

Step  4b.)  The  existing  feasible  flow,  as  illustrated 
in  Figure  8,  is  optimal.  The  associated  optimal  cost  of 
sending  this  flow  through  the  original  network  is  $46. 
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V.  EXTRUSIONS 


A.  SPECIAL  CASE:  Q  1  MINIMUM  .'1^  VALUE 

Consider  a  problem  where  the  specified  amount  of  flow  Q 
to  bo  sent  through  the  network  is  less  than  or  equal  to  the 
value  of  every  arc  in  the  network.  Essentially,  tho 
problem  presented  by  Zangwill  [8]  falls  under  this  catogory, 
since  he  did  not  even  consider  the  existence  of  upper  bounds 
on  arc  flow  capacities. 

Such  a  problem  can  be  solved  by  using  the  algorithm 
presented  in  this  paper.  However,  an  easier  solution  pro¬ 
cedure  to  this  problem  exists.  A  minimum-cost  chain  from 
the  source  to  the  sink  can  be  identified  by  employing  any 
shortest  route  algorithm,  such  as  those  discussed  by 
Dreyfus  [4] . 

Let  d^j  represent  the  length  of  the  directed  arc  from 
nods  i  to  node  j  in  a  shortest  route  algoritlun,  such  *? 
the  one  developed  by  Dijkstra  [3] .  According  to  Dreyfus, 
Dijkstra's  shortest  route  algorithm  is  computationally  the 
most  efficient. 

For  the  above  problem,  let 


ij(Q>, 

oo  , 


if  arc  (i,j)  exists  in  the  network 

(7) 

if  arc  (i,j)  docs  not  exist  in 
the  network 


To  use  Dijkstra's  algorithm,  all  d^j  values  must  be 
greater  than  or  equal  to  zero.  Since  tho  cost  functions 


25 


are  uBuuinocl  to  bo  non-negative  over  the  entire  range  of 

their  respective  foaulblo  arc  flown,  thin  com! i  Lion  in  met. 
(k) 

Also,  let  denote  tho  label  value  of  nodo  i  after 

the  kth  iteration  of  bijkstra'u  algorithm  ban  boon  completed. 
Dijkstra'u  alcjorithin  procouda  uu  follows: 

Stop  0.)  Sot  v|0)  •  0,  V<°>  -  oo  ,  i  «  2,3,... ,n. 

Doclaro  ■  0  as  tho  permanent  label  valuo  of 

node  1. 

Step  1.)  Compute  m  min  |*vj0^  ,  •  Then  find 

V<X>-  min  and  specify  V,  ■  ou  tho  permanent  label 

P  j  P  p 

valuo  of  node  p. 

Stop  2.)  Compute  ^  “  min  VP+l,pj}  '  Than 

find  »  min  V.j2*  and  specify  VtJ  =  Vtj2^  a»  the  permanent 

label  value  of  node  q. 

Step  3.)  Continue  the  implied  iterative  process  until 
node  n  has  a  permanent  label  value. 

At  most,  n-1  iterations  will  be  required  to  label  node  n. 
The  minimum-cost  chain  will  consist  of  the  arcs  (l,p),  (p,q) , 
(q,r),...,  (m,n).  Send  Q  units  of  flow  along  this  chain. 

B.  NON-ZERO  LOWER  BOUNDS 

Tho  algorithm  presented  in  this  paper  can  bo  modified 
to  handle  problems  in  which  non-zero  lower  bounds  are 
associated  with  the  arcs  of  the  network.  Let  L^j  be  the 
lower  bound  on  flow  capacity  for  an  arc  (i,j).  Then  for 
X^j  to  be  a  feasible  arc  flow, 
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'ij  -  xi) 


i  M 


ij 


(8) 


for  each  arc  (i,j)  in  the  network. 

Constraint  (6)  in  the  cost-minimization  problem  is 
replaced  hy  (0).  Once  again,  the  algoritlim  begins  by 
finding  any  flow  which  satisfies  the  constraints.  Thu 
only  change  in  the  algorithm  is  that  a  length  sequence 
network  is  constructed  no  follows: 

1. )  If  arc  (i,j)  of  tho  original  network  is  at  its 
lower  bound,  draw  a  forward  arc  from  node;  i  to  mnlo  j  in 
tho  length  sequences  network.  Label  this  forward  arc  with 

,<D  -  <Mi  i-'-i  j> 

the  length  sequence  ^ ij  *  wh°ro 

(k)  (L. . +k) 

I  a  -  cuiJ  • 

2. )  If  arc  (i,j)  of  the  original  network  is  neither 
at  its  lower  bound  nor  saturated,  draw  a  forward  arc  from 
node  i  to  nodo  j  and  draw  a  backward  arc  from  node  j  to 
node  i  in  the  longth  sequence  network.  Label  the  forward 

(1)  <Mirxil> 

arc  with  the  length  sequence  I  ^  J  ,  whore 

(k)  (XA1+k) 

JPAj  ■  J  .  Label  the  backward  arc  with  the  length 

(1)  ,(xii"Lii)  l(k>  (X,  +l-k) 

sequence  ,...,  | ^  J  where  "ji  *  J 

3. )  If  are  (i,j)  of  the  original  network  in  saturated, 
draw  a  backward  arc  from  nodo  j  to  node  i  in  the  length 

soquonce  network.  Label  this  backward  arc  with  the  length 

(1)  <M  -L  i 

sequence  I4l  j  13  ,  whore 
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(k) 


ji  "  ’Cij 


ji 

(M^+l-k) 


ji 
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vi.  ah i : as  i'ui<  j’manu^nTuov 


a.  identification  ok  negative  cycu.n 


The  only  drawback  with  tho  algorithm  preuenUd  in  thin 
paper  in  that  it  includes  no  procedure  lor  idrnt  i  fyiny 
cycles  of  negative  length  in  tho  length  sequence  notvork. 
Since  tho  identification  of  nuch  cycli  c  ycncratcr.  the 
iterative  pi  ocean  of  the  algorithm,  tho  discovery  of  an 
organized  procedure  which  identifies  such  cycles  would  be 
most  useful. 


Yon  17)  has  devised  an  efficient  shortcut  route  al¬ 
gorithm  which  makes  forward  and  backward  scanning  pannes 
on  a  matrix  coni>oncd  of  values  as  defined  above.  Au  a 
by-product  of  his  algorithm,  simple  cycles  of  negative 
lenyth  can  easily  Ihj  identified.  The  attempts  by  the 
author  to  identify  multiple  and  compound  cycles  by  modifying 
Yen's  matrix  were  not  successful. 


B.  GENERAL  COST  FUNCTIONS 

The  ulgoritl>m  presented  in  this  paper  can  also  be  used 
to  solve  problems  in  which  the  cost  functions  arc  linear 
or  convex.  Due  to  tho  non-decreasing  naturo  of  the  related 
coot  increments  and  length  sequences,  it  would  only  bo 
nocoesary  to  find  simplo  cycles  of  negative  length.  The 
restriction  against  using  both  members  of  a  complement 
could  be  waived. 
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It  is  the  author's  contention  that  the  algorithm  could 
be  modified  to  solve  problems  in  which  the  cost  functions 
are  neither  convex  nor  concave.  Further  research  in  this 
area  might  prove  to  be  very  interesting. 
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VII.  SUMMARY 


A  procedure  has  been  developed  and  presented  for 
optimally  sending  any  specified  amount  of  flow  through  a 
network  in  which  each  arc  has  an  upper  bound  on  capacity 
and  a  concave  cost  function. 

The  solution  technique  involves  the  construction  of  a 
length  sequence  network  corresponding  to  a  feasible  solution 
of  the  problem.  The  identification  of  negative  cycles 
in  the  length  sequence  network  results  in  a  recirculation 
of  flow  in  the  original  network  and  yields  a  reduction  in 
the  overall  cost.  Recirculation  of  flow  in  the  original 
network  creates  another  length  sequence  network.  The 
iterative  process  continues  until  no  further  negative  cycles 
can  be  identified.  When  this  occurs,  the  existing  feasible 
flow  in  the  original  network  is  optimal. 

Inputs  needed  to  solve  the  problem  are  the  upper  bound 
on  capacity  and  the  cost  function  associated  with  each  arc 
in  the  network. 

Modifications  to  the  solution  procedure  are  presented 
so  that  network  problems  having  the  additional  characteristic 
of  non-zero  lower  bounds  may  be  solved.  A  very  simple 
substitute  method  is  also  presented  to  deal  with  a  special 
case  of  the  problem:  Q  1  minimum  M. .  value  in  the  network. 


30 


BIBLIOGRAPHY 


1.  Dantzig,  G.  3.,  Linear  Programming  and  Extensions, 

p.  545,  Princeton  University  Press,  1963. 

2.  Dantzig,  G.  B.,  Ford,  L.  R. ,  and  Fulkerson,  D.  R. , 

A  Primal  Dual  Algorithm,  The  RAND  Corporation, 

Paper  P-778,  5  December  1955. 

3.  Dijkstra,  E.  V7.  ,  "A  Note  on  Two  Problems  in  Connection 

with  Graphs,"  Numerische  Mathematik,  Vol.  1, 
pp.  269-271,  1959 . 

4.  Dreyfus,  S.  E. ,  "An  Appraisal  of  Some  Shor tost-Path 

Algorithms,"  Operations  Research,  Vol.  17,  No.  3, 

May- June  1969’. 

5.  Ford,  L.  R.  and  Fulkerson,  D.  R. ,  Flows  in  Networks, 

Princeton  University  Press,  1963. 

6.  Hu,  T.  C.,  "Minimum-Cost  Flows  in  Convex-Cost.  Networks," 

Naval  Research  Logistics  Quarterly,  Vol.  13,  No.  1, 
March  1966. 

7.  Yen,  J.  J.,  Some  Algorithms  for  Finding  the  Shortest 

Routes  Through  the  General  Networks ,  paper  presented 
at  the  SIAM  2nd  International  Conference  on  Computing 
Methods  in  Optimization  Problems,  San  Remo,  Italy, 

9-13  September  1968. 

8.  Zangwill,  W.  I.,  "Minimum  Concave  Cost  Flows  in  Certain 

Networks,"  Management  Science,  Vol.  14,  No.  7, 

March  1968. 


31 


